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Abstract 

The spontaneous baryogenesis mechanism by Cohen, Kaplan and Nelson is reconsidered taking 
into account the transport of particles inside the electroweak bubble walls. Using linear response 
theory, we calculate the modifications on the thermal averages of the charges of the system 
due to the presence of a space time dependent 'charge potential' for a quantum number not 
orthogonal to baryon number. The local equilibrium configuration is discussed, showing that, 
as a consequence of non zero densities for conserved charges, the B + L density is driven to non 
zero values by sphaleronic processes. 

Solving a rate equation for the baryon number generation, we obtain an expression for 
the final baryon asymmetry of the Universe containing the relevant parameters of the bubble 
wall, i.e. its velocity, its width, and the width of the region in which sphalerons are active. 
Compared to previous estimates in which transport effects were not taken into account, we find 
an enhancement of nearly three orders of magnitude in the baryon asymmetry. 

Finally, the role of QCD sphalerons in cooperation with transport effects is analyzed, show- 
ing that the net result depends crucially on the particle species which enter into the charge 
potential. 
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1. Introduction 



The possibility of a baryogenesis at the electroweak scale is a very popular but controversial 
topic. Despite the large number of related publications none of the key aspects of this 
subject can be considered to be on firm ground. First, it is well known that the requirement 
that the anomalous 'sphaleronic' processes which violate baryon number (B) go out of equilib- 
rium soon after the transition translates into a lower bound on the ratio between the vacuum 
expectation value (VEV) of the Higgs field, v{T), and the critical temperature of the transition, 
Tc, f (Tc)/Tc> 1. In the standard model, improved perturbative evaluations of this ratio give 
a value which is badly less than unity for values of the mass of the Higgs scalar compatible with 
LEP results. The situation in the minimal supersymmetric standard model is slightly better, 
however the allowed region of the parameter space is very small, and is likely to be excluded 
by future LEP data On the other hand, the perturbative expansion cannot be trusted 
any more for values of the Higgs masses of the order of the W boson mass or larger 0], and 
preliminary results based on non perturbative methods (lattice simulations |Q, 1/e expansion 
P, effective action ^) give indications of strong differences of the results with respect to those 
obtained perturbatively. Clearly, much work is still needed on this issue. 

Another aspect of the problem is CP violation. This has been the subject of a recent debate 
in the literature about the need of further complex phases in the theory besides the one in the 
Cabibbo-Kobayashi-Maskawa mixing matrix (see. |^). In models with more than one Higgs 
doublet, like the minimal supersymmetric standard model, further sources of CP violation can 
emerge naturally from the Higgs sector. In particular, the possibility of a spontaneous CP 
violation at finite temperature has been emphasized 10, 11 1. This effect could give enough 
contribution for the baryogenesis and at the same time satisfy the upper bounds on CP violation 
coming from the electric dipole moment of the neutron. 

Even assuming that the phase transition is strong and that CP violation is enough, we 
must face the other key issue, namely what is the mechanism responsible for the generation 
of baryons during the phase transition. The most significant departure from thermodynamic 
equilibrium takes place at the passage of the walls of the expanding bubbles which convert 
the unbroken into the broken phase. According to the size and speed of the bubble walls, two 
different mechanisms are thought to be dominant. In the case of "thin" (width ~ ^/T) walls, 
typical of a very strong phase transition, the creation of baryons occurs via the asymmetric 
(in baryon number) reflection of quarks from the bubble wall, which biases the sphaleronic 
transitions in the region in front of the expanding bubble [|T^ . If the walls are "thick" (width~ 
(10 — 100)/T) then the relevant mechanism takes place inside the bubble walls rather than 
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in front of them. In this case we can make a distinction between 'fast' processes (mediated 
by gauge, flavour diagonal, interactions and by top Yukawa interactions) and 'slow' processes 
(mediated by Cabibbo suppressed gauge interactions and by light quarks Yukawa interactions). 
The former are able to follow adiabatically the changing of the Higgs VEV inside the bubble 
wall, while, in first approximation, the latter are frozen during the passage of the wall. If CP 
violation, explicit or spontaneous, is present in the scalar sector then a space-time dependent 
phase for the Higgs VEVs is turned on inside the wall. The time derivative of this phase 
couples with the density of a quantum number non orthogonal to baryon number (e.g. fermion 
hypercharge density) ^ and then can be seen as an effective chemical potential, named "charge 
potential" , which has the effect of biasing the rates of the sphaleronic processes, creating an 
asymmetry proportional to -(9, where i} is the phase of the VEVs. 

This "adiabatic scenario" , originally due to Cohen, Kaplan, and Nelson (CKN) has been 
recently reconsidered by different authors in different but related aspects. First, Giudice and 
Shaposhnikov have shown the dramatic effect of non perturbative, chirality breaking, transitions 
induced by the so called "QCD sphalerons" . If these processes were active inside the bubble 
walls, then the equilibrium value for baryon number in the adiabatic approximation would be 
proportional to that for the conserved quantum number B — L (L is the lepton number), up 
to mass effects suppressed by ~ {mtop{T) / tiT)"^ . Then, imposing the constraint {B — L) = 
(here (■ ■ ■) represents the thermal average) we obtain zero baryon number (up to mass effects). 



Dine and Thomas have considered the two Higgs doublets model in which the same 
doublet couples both to up and down quarks, the same model considered in the original work 
by CKN. These authors have pointed out that ^ couples also to the Higgs density, so that the 
induced charge potential is for total hypercharge rather than for fermion hypercharge. As long 
as effects proportional to the temperature dependent VEV v(T) are neglected, hypercharge 
is a exactly conserved quantum number and then, again imposing the constraint that all the 
conserved charges have zero thermal averages, no baryon asymmetry can be generated. So, we 
again find a mtopiT)"^ /T"^ suppression factor. 

Finally, Joyce, Prokopec and Turok (JPT) have emphasized the very important point that 
the response of the plasma to the charge potential induced by d is not simply that of a system 
of fixed charges, because transport phenomena may play a crucial role |17]. When a space time 
dependent charge potential is turned on at a certain point, hypercharged particles are displaced 
from the surrounding regions, so that even the thermal averages of conserved quantum numbers 
become locally non vanishing. As a consequence, the equilibrium properties of the system have 
to be reconsidered taking into account the local 'violation' of the conserved quantum numbers. 



^The phase also couples to Chern-Simons number but this coupling induces an effect which is suppressed by 
Wtop(r)^/T^ with respect to the one which we are presently discussing ll^, so we will neglect it. 
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In this paper we analyze the adiabatic scenario using hnear response theory |jTB|, in order 
to take transport effects into account. We assume that a spacetime dependent charge potential 
for fermion hypercharge is generated inside the bubble wall, without discussing its origin, and 
investigate its effect on the thermal averages of the various quantum numbers of the system. 
We find that transport phenomena are really crucial, but we disagree with JPT's conclusion 
that as a consequence of the local 'violation' of global quantum numbers there is no biasing 
of the sphaleronic processes. Actually, in the adiabatic approximation the local equilibrium 
configuration of the system is determined by the thermal averages of the charges conserved by 
all the fast interactions. The effect of transport phenomena is to induce space time dependent 
non zero values for these averages. We calculate these averages using linear response theory and 
then determine the local equilibrium configuration, showing that it corresponds to (B + L) ^ 0. 
As we will discuss, JPT's result corresponds to freezing out any interaction inside the bubble 
wall, which is in contradiction with the adiabatic hypothesis. Then we write down a rate 
equation in order to take into account the slowness of the sphaleron transitions and obtain 
an expression for the final baryon asymmetry explicitly containing the parameters describing 
the bubble wall, such as its velocity, v^^ its width, and the width of the region in which the 
sphalerons are active. 

The inclusion of transport phenomena also sheds a new light on the strong sphaleron ef- 
fects and on the effect of a charge potential for total rather than fermionic hypercharge. The 
dramatic suppressions found by Giudice and Shaposhnikov and by Dine and Thomas respec- 
tively, are both a consequence of taking zero averages for conserved quantum numbers. Since 
these averages are no more locally zero we will find a non zero {B + L) ^ 0, even in the case 
in which the charge potential is for total rather than for fermion hypercharge. In the case 
of QCD sphalerons we find that the final result depends in a crucial way on the form of the 
charge potential which is considered. For example, if all the left handed fermions plus the right 
handed quarks contributed to the charge potential according to their hypercharge, then no bias 
of sphaleron processes would be obtained. In this case, we would find a non zero value for B+L 
inside the bubble wall but no baryon asymmetry far from it inside the broken phase. On the 
other hand, if only right and left handed top quarks participate to the charge potential, then a 
final asymmetry is found and QCD sphalerons are harmless. 

The plan of the paper is as follows: in section 2 we will develop a chemical potential 
analysis for the equilibrium properties of the plasma inside the bubble wall in the adiabatic 
approximation, and we will write down the rate equation fro the production oi B + L due to 
sphaleron transitions; in sect. 3 we will introduce our application of linear response analysis 
to the calculation of the variations of the thermal averages induced by d. In particular we will 
show that the fundamental quantity to evaluate is the retarded two points Green's function for 
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fermion currents. In section 4 we will solve the rate equation, finding an expression for the final 
baryon asymmetry in terms of the various bubble wall parameters. In this context, we will also 
discuss the screening effects on the electric charge. In sect. 5 we will discuss the role of QCD 
sphalerons in cooperation with transport phenomena, and finally we will discuss our results in 
section 6. 



2. Local equilibrium inside the wall 



For definiteness, we will work in the two Higgs doublet model in which one doublet couples 
to up quarks and the other one to down quarks. The phase transition is assumed to be strong 
enough so that sphaleron processes freeze out somewhere inside the bubble wall (see the end 
of sect. 4 and ref. [1^ for a discussion about this point). 

The relevant timescale for baryogenesis is given by the passage of the bubble wall, which 
takes Atu) = Az/vy, ~ (200)/T where Az is the wall thickness. During this time the 
phase of the Higgs VEV's changes of an amount A-t}. Thus, we can discriminate between fast 
processes, which have a rate >1/At^t, and then can equilibrate adiabatically with and slow 
interactions, which feel that is changing only when the bubble has already passed by and 
sphalerons are no more active. Next, we introduce a chemical potential for any particle which 
takes part to fast processes, and then reduce the number of linearly independent chemical 
potentials by solving the corresponding system of equations, in a way completely analogous to 
that followed for example in refs. |2^, the main difference here being that light quark Yukawa 
interactions and Cabibbo suppressed gauge interactions are out of equilibrium. Finally, we 
can express the abundances of any particle in equilibrium in terms of the remaining linear 
independent chemical potentials, corresponding to the conserved charges of the system. 

Since strong interactions are in equilibrium inside the bubble wall, and since the current 
coupled to is color singlet, we can chose the same chemical potential for quarks of the same 
flavour but different color, and set to zero the chemical potential for gluons. Moreover, since 
inside the bubble wall SU(2) x f/(l) is broken, the chemical potential for the neutral Higgs 
scalars vanishes^. 

The other fast processes, and the corresponding chemical potential equations are: 
i) top Yukawa: 



^This is true if chirality flip interactions, or processes like Z Z*h, are sufficiently fast; since the corre- 
sponding rates depend on the Higgs VEV, they will be suppressed by factors of {v{T)/TY with respect for 
example to the rate for Ml ^ tug. This has led the authors of refs.[^ to consider the system in the 
unbroken phase. Anyway, this choice does not lead to dramatic changes to the conclusion of this paper. 
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ii) SU(2) flavour diagonal: 



UL*^dL + W = ^^d\ + ^^W+) / • _ T 9 ox /ON 

rrO 7T+ , TI/" ( \ I* " J-, ^5 <jj l^J 

^H- + W+ [fiH- = -l^w+) 

Neutral current gauge interactions are also in equilibrium, so we have zero chemical potential 
for the photon and the Z boson. 

Imposing the above constraints, we can reduce the number of independent chemical po- 
tentials to four, Hw+i l^tLi I^UL = 1/2 Z]i=i f^eL = 1/3 Z^Li /^e*^- These quantities 
correspond to the four linearly independent conserved charges of the system. Choosing the 
basis Q', {B - L)\ {B + L)', and BP' = B'^ - 1/2{B[ + B'^), where the primes indicate that 
only particles in equilibrium contribute to the various charges, and introducing the respective 
chemical potentials, we can go to the new basis using the relations 




(3) 



= 3/if^ — 2fiuL 

If sphaleron transitions were fast, then we could eliminate a further chemical potential 
through the constraint 

3 3 3 

2E/^«l+E/^4+E/^ei =0- (4) 

i=l i=l i=l 

In this case, the value of {B + L)' would be determined by that of the other three charges 
according to the relation 

3 7 19 

(B + L)' =—Q' + —BP' -—(B-L)' (5) 

where we have indicated charge densities by the corresponding charge symbols. We have ne- 
glected mass effects, which means that the excess of particle over antiparticle density is related 
to the chemical potentials according to the relation |20| n+ — n_ = agT^/6{fi/T), where a = 1 



for fermions and a = 2 for bosons, and g is the number of spin and color degrees of freedom. 



The above result should not come as a surprise, since we already know from ref. that a 
non zero value for B — L gives rise to a non zero i? + L at equilibrium. Stated in other words, 
sphaleron transitions erase the baryon asymmetry only if any conserved charge of the system 
has vanishing thermal average, otherwise the equilibrium point lies at {B + L)eq 0. 

Actually, sphaleron rates are too small to allow {B + L)' to reach its equilibrium value (^), 
Tgp ~ (a^T)^^ ^ ^tw: so equilibrium thermodynamics cannot be used to describe baryon 
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number generation inside the bubble wall. Following refs. |T1| we shall make use of the rate 
equation 

^ (u \ TV (a\ 

d^B + ^^SP = -^g^^^ (6) 

where Tsp = k{awT)^ exp{—(f)/gwT) is the rate of the sphaleron transitions when the value 
of the Higgs field is (/c ^ (0.1 — 1) from numerical simulations [^), F is the free energy of 
the system, and the derivative with respect to {B + L)' must be taken keeping Q', {B — L)', 
and BP', constant. The meaning of eq. (^) is straightforward. Sphaleron transitions (which 
change {B + L)' but conserve Q', {B — L)', and BP') will be turned on only if they allow the 
total free energy of the system to get closer to its minimum, i.e. equilibrium, value. At high 
temperature (/i^ <^ T) the free energy of the system is given by 



F =^ 

^ 12 



Using ([l|) , d^) and (^ to express the chemical potentials in terms of the four conserved charges 
in (I) we obtain the free energy as a function of the density of {B + L)' = (ib+lT^/G, 



F[{B + L)'] = 0.46 



{B + L)' -{B + L)'^Q ' 



2^2 



+ constant terms 



where the "constant terms" depend on Q' , {B — L)', and BP' but not on {B + L)', and {B + L)'^q 
is given by (|^). The total amount of {B + L)' present in a certain point at a certain time is 
made up by two contributions: {B + L)'gp, generated by sphaleron transitions, and {B + L)'rp^, 
which is not generated but is transported from nearby regions in response to the perturbation 
introduced by 'd. So, eq. (||) now takes the form 

j^{B + L)'sp = -0.92^ [{B + L)'sp + {B + 1)'^^ - (B + (9) 

with the initial condition {B + L)'gp = before i) is turned on. 

Let us summarize our discussion up to this point. Consider an observer in the plasma 
reference frame during the phase transition. When a bubble wall passes by the observer, he 
measures a space time dependent charge potential for, say, fermionic hypercharge, which induces 
transport phenomena and then local asymmetries in particle numbers. The Q', [B — L)', and 
BP' components of these asymmetries remain unaffected by fast interactions, while the other 
components are reprocessed as to obtain their equilibrium values. In the case of {B + L)' the 
reprocessing is slow, so we have to use the rate equation in (^ to describe it. The generation of 
[B + L)'gp will go on until either Tsp goes to zero or the local equilibrium value is reached i.e. 
{B + L)'sp + {B + L)'rpj^ = {B + L)'p,Q. After the passage of the wall, {B + and {B + L)!^^ 
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go rapidly to zero since ^ vanishes, and so the final asymmetry is given by the {B + Ly^p 
generated until that time. 

As we can see, the crucial question is now to calculate the induced values for Q', {B — L)', 
BP', and {B + L)'j,f^. We will do that in the next section by using linear response analysis 



18, pa]. 



3. Linear response analysis 



In this paragraph we want to discuss the effect on the thermal averages of the term induced 
in the lagrangian when d is active, which we assume to have the form 

C^C + §J%, (10) 

with 

= E/yW- (11) 

where means that the sum extends on particles in equilibrium with only. The standard 
procedure fl^ is to consider as an effective chemical potential, so that particle abundances 
at equilibrium are given by 

Pi = liPQ + (b- l)iHB-L + bpifiBP + yf^, (12) 

and then imposing 

{Q') = {{B - L)') = {BP') = (13) 

so that any particle abundance can be expressed in terms of only. The point is that taking 
transport phenomena into account, the above thermal averages are not zero, but depend on -(9 
themselves. So we must first calculate their values and then use them to determine the chemical 
potentials. Our starting point is the generating functional 

Z[Jo;^]= f V(l)exp(i f dr f £x[C + ^J^^ + JoO + ' sources']] (14) 

where T><^ is the integration measure, O is the operator of which we want to calculate the 
thermal average and Jq the corresponding source. C is any path in the complex r plane 
connecting the point rj„ to Tout = '^in — iP {P = '^/T) such that the imaginary part of r is never 
increasing on the path p4|. P(A)BC means that periodic (antiperiodic) boundary conditions 
must be imposed on bosonic (fermionic) fields on the path. 

The thermal average of the operator 0(r, x) in presence of is obtained in the usual way 

1 SZ[Jo;^] 



(0(t,,x)) 



i 6Jo{tx,x) 



(15) 
Jo=o 
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'Tin tx 



nir) 



Tout 



7^1(r) 



-i/3 



Figure 1 : The path C corresponding to imaginary time formahsm of thermal field theory. 



where all the field sources are set to zero. Now we make a functional expansion of (|T^) in 
and truncate it at the linear term, 
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i 5Jo{t^,x) 



Jc Jv ^ ^ Sd{r',y) 



i9=0 



(16) 



Jo=o 



{Oitx,x))^__ 



■-Jo=o 



where Tc is the ordering along the path C. 

Different choices for the 'time' contour C lead to different formulations of thermal field 
theory. One possibility is to take the vertical line connecting tx to tx — so that r' — tx is 
pure imaginary on any point of the path. 

This choice corresponds to the imaginary time formalism (ITF) of thermal field theory 
and in this case we have to evaluate the euclidean two point thermal Green's function 



ZE),om{z 



This can be done in Matsubara formalism, where FejTiman 



rules are straightforwardly obtained ||2^. But, in this case, we would have to calculate "d for 
complex times, whereas we are interested in its values at real times, during the passage of the 
wall. So, in order to get a more direct physical interpretation of what we are calculating, we 
must turn to real time formalism. This corresponds to choose the path in Fig. 2, and then 
letting Tin going to — oo, and tp to +oo Now, it is possible to see that the contributions 
to the integral coming from r' on C3 vanishes in the above limit |2^, so we are left with the 
contributions from Ci and C2 only. The Tq ordering now allows us to rewrite the integral in 
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Tin tx 



nir)' 



Co 



Ci tj 



7^ll 



Tout 



-ip 



Figure 2: The path C corresponding to real time formulation of thermal field theory. C2 lies infinitesimally 
beneath the real axis. 



(|T6]) as 



i dty d^y^{ty,y) {TcJYpity,y)0{t^,x))^^j^^f^ 

J Ci®C2 •'V 

= i f dty f d^y^{ty,y) {\jY^{ty,y),0{t^,x)])^^j 



(17) 



Defining as usual the retarded Green's function 



iDovpit^,x;ty,y) = { 0{t^,x), Jy^ity^y) )Q{t^ - ty) 



where 0(x) is the step function, we arrive at the result 



t-oo 



(O(t.,x))^.o = (O(t,,x))^^o+ / dty d-'y^{ty,y)D^y^{t,,x-ty,y). (19) 

J —00 JV 

The operators we are interested in are fermionic charge densities {Q\ {B — L)', {B + L)', BP') 
of the form Q'^ = Y^^qfJi- Inserting it in (|^), and using definition (|ll]) we get 

/+00 r 
dty / d^y§{ty,y)Df^{t^,x-ty,y) (20) 
-00 JV 

where Df^ is the current-current retarded Green's functions for fermion i and j [i and j are 
flavour and color indices) and we used the fact that {Q'A(tx,x))^^Q = 0. Note that the Green's 
function has to be evaluated for i) = 0, i.e. we must use the unperturbed lagrangian with 
the charge potential turned off and all the chemical potentials equal to zero in the partition 
function. 

The problem of calculating the effect of the charge potential in (|10[) on the thermal averages 
for the particles in equilibrium is then reduced to the evaluation of the retarded Green's func- 
tions which enter in (pOl). As we have discussed, the more natural framework for this calculation 
is real time thermal field theory, in which the physical sense of the various quantities is evident. 
Anyway, we have also seen that we can calculate the euclidean Green's function in the imagi- 
nary time formalism and then continue analytically to real times, thus obtaining the -Df^'s (this 
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relation was established for the first time by Baym and Mermyn [^). In energy- momentum 
space the analytical continuation is accomplished by the substitution 

iujn uj + ie e ^ (21) 

where a;„ = 27mT are Matsubara frequencies and u is the real energy. 



4. Solution of the rate equation 



Since the rate of the sphaleronic transition is suppressed by a^, the asymmetry in {B + L)' 
generated by the sphalerons, {B + Ly^p, is generally much smaller than both {B + LY^q and 
[B + L)^^ (we can check it a posteriori) , so we can approximate the rate equation in (^ by 

+ LYsp ^ 0.92^ [{B + LYeq -{B + L)^ . (22) 
Using the equations (||) and (^Of ) we can determine {B + LYi^^ = (5 + L)^g — (5 + L)^^ as 

/+00 r 
dty / dV^(i„ y) Df^{t., X- ty, y) (23) 
-oo JV 



where 



3 7 , 19 



.80' 20 ' 40 

Integrating eq. ( p2l) in time from — oo to +00 we get the final density of {B + L)', 



/+00 ^ 
rft, — (5 + L)'sp(t,,f) 

271" /•+00 _ / 

= 0.92— y du;rsp(-u;,x) (5 + L)^^(a;,x) 



(24) 



where f sp{—^-,x) and {B + L)^p(uj,x) are the Fourier transformed of r5p(ta;,x) and (B + 
L)LR{tx,x) with respect to time. 

We recall that Tsp is /c(aiy7')'* in the unbroken phase and decreases exponentially fast as the 
Higgs VEV is turned on. In order to solve eq. (|2^ ) analytically we approximate this behaviour 
by a step function. Moreover, we will consider a plane bubble wall moving along the z-axis 
with velocity Vw So, our expression for Tsp will be 

Tspit,, f) ~ r e ft, - ti - ^) e L - + — ) (25) 

with ti —00 and F = k{awT)'^. Of course, more sophisticated approximations for Tsp may 
be used, at the price of solving eq. (E3) numerically. 
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X = 



tl 



I r 



Figure 3: Schematic representation of the behaviour of d and of the rate of the sphaleronic transitions at the 
point X = 0. 

Analogously, we approximate in such a way that it is constant in a region of width Az 
inside the bubble wall, and is zero outside, 

-dity, y) =6 Q{Zy - Vyjty) Q{Vyjty ^ Zy + Az) (26) 

where 9 = Vw A'd/Az. So, if we are at the point x = 0, we observe an interaction of the form 
( p!0D turned on from t = —Az/v^ to t = 0, while the sphalerons are active till t = ^2, as we 
have shown in Fig. 3. 

Putting all together we obtain 

(27) 

du Dj^jiPx = Py = 0,Pz = ,t^)- 

oo io u; Vw 

Note the peculiar relationship between pz and u in the argument of Df^, due to the spacetime 
dependence of 'd(ty,y), see (p6|) . A consideration of the general properties of retarded Green's 
functions ||26| ensures that the imaginary part of is a even function of u, while its imaginary 
part is odd. As a consequence the integral in (^7]) will always give a real result. 

The lowest order contribution to Dj^ comes from the fermion loop in Fig. 4, where the two 
crosses indicate the zero components of the fermion current. We may evaluate the corresponding 
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J? 



Figure 4: Lowest order contribution to D 



euclidean two point Green's function in the ITF and then continue analytically to real energies 
according to eq. (|2TD . Moreover, since the relevant frequencies at which Df- must be evaluated 
are such that uj <\p\ < v^j IS.z ^ we take only the leading terms in the high temperature 

expansion. These are given by pTl 



ni(po = 2n7rT,p) bij 



3 



2\p\ ipo - \p\ 



(28) 



where we have neglected fermion masses. After continuing analytically to real energies and 



fixing the momenta as in (|27|) we get the lowest order contribution to Df^, 



Py = 0,p, = —,UJ) = — 



1 



2 1 — f „ 



TT 



+ i-v.^sign{uj^ 



5ij . 



(29) 



When f ^ — >■ 1 the above expression exhibits a collinear divergence, due to the fact that the 
fermions in the loops are massless in our approximation. This divergence disappears when 
plasma masses for fermions are taken into account. However, for our purposes, since — 0.2 



, the effects of plasma masses for fermions can be neglected [|28 
Due to the constraint p 



ujjvyj the real part of the correlation function in (p9D does not 
depend on cj, while the imaginary part depends on its sign only. This implies that the response 



induced on the plasma through (^9]) has neither spatial nor temporal dispersion, i.e. inserting 
( p9D in (pOD gives rise to an induced thermal average for the charge Qa which in any space-time 
point is proportional to the value of d in that point 



(30) 



In particular also {B + L)j'r and {B + L)eq receive a contribution of this form and disappear as 



soon as 'd is turned off. Inserting it into the rate equation we obtain from (pTj) the contribution 
to the asymmetry {ti — oo) 



A(i? + L)', 



SP 



(31) 
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Figure 5: The contribution to D(j due to photon exchange. 

where 



< 

+ < ^2 < O) (32) 



Az 



We recall that t2 is the time at which sphaleron transitions are turned off, while the charge 
potential induced by is active for —Az/v^, < t < 0. So the asymmetry calculated in this 
approximation for Dj^j grows linearly with t2 until ^2 = (for t2 < —Az/v^ the asymmetry is 
obviously zero since there is no overlap between sphalerons and The result for ^2 > is an 
artifact of our approximation {B + L)'gp ^ (S + L)'j^^, which is no more appropriate in this 
case. Actually, from (^) we know that when goes to zero, as is the case for t > 0, {B + L)'^q 
and {B + L)'j.^ vanish too, and the rate equation (j^) becomes 

j^{B + lysp = -0.92^ [{B + L)sp\ {t > 0) (33) 

so that the asymmetry produced before decreases exponentially from t = to t = ^2 with rate 
r. However, due to the smallness of F, and to the fact that t2 cannot be much larger than 
0{Az/vyj)^ we can safely neglect this decreasing and take the result in (p^). 

Next, we consider the contribution to Df^ due to gauge bosons exchange. Since we are 
in the broken phase, and since the perturbation ( pUj ) induced by d is colorless, we will take 
into account only photons, which contribute through the graph in Fig. 5. The blob in the 
middle represents the sum of all possible insertion of fermion loops. In calculating the blob, 
we have to include not only the fermions which enter in Q\ {B — L)', and BP', i.e. the ones 
with fast flavour, or chirality, changing interactions, but we must instead take into account the 
contributions of all the charged fermions of the theory. This is because QED interactions are 
fast and then, for instance, pair production is in equilibrium also for right handed light quarks. 
The 'QED' contribution of Fig. 5 gives 

I5--(..P1 ^ e%,,Un.,pl^^£^gL__ (34) 
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where eqi is the electric charge of the fermion i. V^^{uj,p) is the tree level (0,0) component of 
the photon propagator in the Coulomb gauge 

= zlpM^ _ A^uV , (35) 
where = uj"^ — = (1, 0, 0, 0) identifies the plasma reference frame, while 

poo ^ pOi ^ pip 

P^' = S,,-pY/\p\' 

As we have discussed, the sum in the denominator of (0) must be extended over all the charged 
quarks and leptons. Setting = u/vw we get 

D^-^^^{u,p) = -e\g, vl , y (36) 



Note that unlike the 'direct' contribution ([29| ) the 'QED' one is not flavour (or color) diag- 



onal, so that even particles which do not enter into the expression for the charge potential (|TOD 
get a non zero thermal average depending on Moreover, this contribution has a genuine oo 
dependence. Two points retarded Green's function are analytic in the upper half of the complex 
u plane. D^j''^^^{u,pz = uj/vw) may then have poles of the form ZJ = Up — i'jp, with > 0. 
In order to determine them we have to solve the following equations 



CltMuJ) (37) 
27p 



Up 



where C = vl^J^ki^Qk) - Since the RHS of the first of eqs. (|37| ) is negative, and v^, — 0.2 < 1, 
we can approximate the solutions by 



7p ^ [Cn\Yii{u)f'' c^vJ^(Y.{qk)' 



1/2 



(3^ 



V3 



k 



where 7p > as it should be. 

Inserting (^) in eq. ( pT]) we obtain the 'QED' contribution to the asymmetry 

A(5 + L)'^p«^^ ^ 0.92^|£r^E:,c., J^^^(t2) (39) 
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where, now, 



1 

3 



2 1 - f „ 



mi 



r 



X 



t2 H h COS U2 H 



Az 



7p 



cos ujpt2 — e 



COS UJp[t2 + 



Az 



it2 < - — 



< t2 < 

1. Vw 



{t2 > 0), 



(40) 

The same considerations about the case ^2 > made after eq. ( p2D apply also here. We can 
notice that photon exchange gives two different types of contributions. The first one has the 
same behaviour of that in (|32|), i.e. a linearly increasing asymmetry from t2 = —Azjv^ to 
^2 = 0. On the other hand, the second term exhibits a well known feature of plasma physics, 
namely, plasma damped oscillations induced by an external perturbation [Q. The damping 
rate here is given by 7p. In the case —Az/v^ < ^2 < 0, we see that the oscillating term 
dominates over the linear one only for t2 — > — A^/t>^, and is rapidly damped as ^2 —^0, since 
exp(— 7pA2;/f^) ^ exp(— TA2;) ^ exp— (40). When t2 > the amplitude of the oscillation is 
always suppressed at least by a factor 10^^ ^ 10^^ with respect to the linear one, and then 
we can conclude that the effect of the oscillating term is negligible unless ^2 is very near to 
-Az/vw. 

An interesting feature of our results (^) and (1^) can be appreciated if we calculate, by 
means of eq. (0), the electric charge Q' induced by the phase ■(9, taking into account both the 
direct contribution (^) and the 'QED' one (|3^) to Df,. It is easy to see that it is given by 



{Q') oc J2i UiQi (Sfc Qk " Sfc Qk) ^ ''linear contribution^' 
+ ^^damped contribution" , 



(41) 



then, when every fermion is in equilibrium, the linear contribution to the thermal average of Q' 
vanishes, and we are left with the damped one. The reason is that in this case Q' coincides with 
the total fermion electric charge, and this is perfectly screened as in the usual QED plasma. 
Since in the real situation not all the fermions are in equilibrium, the linear contribution to (^T]) 
does not cancel. Anyway, this considerations are valid for electric charge only, while the other 
interesting charges, (B + L)', {B — L)', and, BP', would have non zero linear contributions even 
if all the fermions were in equilibrium. 
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Putting all together, and neglecting the damped contribution, we find the final asymmetry 
in {B + L)' 

where we have assumed that the sphalerons turn off when the phase is still active {—Az/vw < 
t2 < 0). Recalling that 6 = Ai^v^/ Az, and assuming that sphalerons cease to be active after a 
time interval ^2 + ^z/vyu = fAz/v^ from the turning on of ^ we get 

A{B + L)'sp - -2.3 ■ IQ^kT^aiA^f (43) 



where we have taken the reasonable value — 0.2 ||T9|. The above value is enhanced by nearly 
three orders of magnitude with respect to the original estimates by CKN [14] where transport 
phenomena were not taken into account. 

The predicted baryon asymmetry of the Universe then comes out to be 

^ ~ -IQ-^kAdf. (44) 



k is estimated in the range 0.1 -r- 1 from numerical simulations |22|, while the value of / is still 



an open question. Following Dine and Thomas |16] we chose 

/ ^ ^ ~ 5 ■ 10-2. (45) 
9 

The observed baryon asymmetry, ps/S = (4 -i- 7) ■ 10"^^, can then be reproduced for Ai) ~ 
10^2^ 10^'^, values which can be obtained either by explicit CP violation or by spontaneous CP 
violation at finite temperature ^ without entering in confiict with the experimental bounds 
on the electric dipole moment of the neutron. 

The result (|4^) was obtained considering a charge potential of the form (|T0|) where the sum 
extends on all the left handed fermions plus the right handed quark. Considering the more 
physical situation in which only the top quarks (left and right handed) feel the effect of i) the 
coefficient 37/240 in (|4^) should be changed into 9/32, thus leading to an enhancement of a 
factor 1.8. 



5. The effect of QCD sphalerons. 

It is well known that the axial vector current of QCD has a triangle anomaly, therefore one 
can expect axial charge violation due to topological transitions analogous to the sphaleronic 
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transitions of the electroweak theory. The rate of these processes at high temperature may be 



estimated as 29, 15 



Tstrong = ^ f— ^^T^P = ^M«sT)^ (46) 

3 \aw/ 3 

where ag is the strong fine structure constant, leading to a characteristic time of order 

^strong = ^^^J^- (47) 



Since /c ~ 0.1 -f- 1 p2|, we see that Tstrong is comparable to the time of passage of the bubble 
wall, and might also be smaller. 

Recently, Giudice and Shaposhnikov have analyzed the effect of these 'QCD sphalerons' 
on the adiabatic baryogenesis scenario. They showed that, as long as these transitions are in 
equilibrium and fermion masses are neglected, no baryon asymmetry can be generated. Thus, 
the final result will be suppressed by a factor ~ {mtop{T) / tiTY . In this paragraph we will 
reconsider the issue taking transport phenomena into account. 

The effect of QCD sphalerons may be represented by the operator 

YiUinLuidLd^), (48) 

where i is the generation index. Assuming that these processes are in equilibrium, we get the 
following chemical potentials equation 

3 

^i.^^ul - /^«»^ + l^di - f^d^) = 0. (49) 

i=l 

Eq. ( ^91 ) contains the chemical potentials for all the quarks, and imposes that the total right- 
handed baryon number is equal to the total left-handed one. Using eqs. (|I|) and (0) we can 
rewrite it as 

4/i„^ + jjt^ - jjb^ - 2^dn - 2/i„^ - 3/ivi/+ = 0, (50) 

where ^iuL r = 1/2 Z]^=i /^u^ ^) and /i^^ = 1/2 X]i=i yU^^- One of the three new chemical poten- 
tials, yUfe^, /irf^, and /iu^, can be eliminated using eq. (pOD, while the remaining two correspond 
to two more conserved charges that must be taken into account besides Q' , BP', and {B — L)' 
(now the primes mean that the summation has to be performed on right handed quarks too, 
but not on right handed leptons). We can choose 

^ = E^k-^E«k, (51) 

i=l i=l 

Y ^ hL + tL + ta + \jZu'n-jl{^i + ^i). (52) 

i=l i=l 
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corresponding respectively to A3 and A2 in the notation of ref. |T^. Following the usual 
procedure, we can now express the abundance of any particle number at equilibrium as a linear 
combination of Q', {B — L)', BP', X, and Y. For {B + L)' we obtain the result 

(i? + L)^Q = -^(5-L)', (53) 

to be compared to eq. (|^), which we obtained considering QCD sphalerons out of equilibrium. 
Thus, the equilibrium value for {B + L)' depends only on the the density of [B — L)', in 
agreement with what was obtained in ref. [00. If transport phenomena were not present, as it 
was assumed in ref. [|15|, we could set {B — L)' to zero and then conclude that QCD sphalerons 
allow no non vanishing [B + L)' density, at equilibrium and in the massless approximation. On 
the other hand, including transport effects, we can easily calculate the {B — L)' density induced 
by 1? using eq. (^), and then, through (0), the equilibrium value {B + L)'^q, which, unlike in 



ref. |15|, comes out to be non vanishing inside the bubble wall. However this is not sufficient to 
conclude that we will have a non zero final baryon asymmetry when the bubble wall has passed 
by. As we discussed in Sect. 2., the generation of baryons inside the bubble walls is described 
by eq. (^, with the initial condition {B + Ly^p = 0. Then we must calculate {B + L)'j.i^, i.e. 
the contribution to {B + L)' due to transport. If all the particles in equilibrium participated 
to the charge potential then we would find 

{B + ly^Rit^, x) = {B + LYEQiU, x) (54) 

so that the system would always be on the minimum of the free energy inside the bubble wall 
and there would be no bias of the (electroweak) sphaleronic transitions. As a consequence, 
{B + L)sp would remain zero and no asymmetry would survive after the passage of the wall 
up to fermion mass effects, in agreement with what was find in ref. ||15 . 



On the other hand, including only left and right handed top quarks into the charge potential, 
eq. (0) is no more satisfied and a non zero result for the final baryon asymmetry is recovered. 
In this case, the factor 37/240 in eq. (|^) should be changed to 25/72. 

6. Conclusions 



In this paper we have analyzed the impact of transport phenomena on the so called 'spon- 
taneous baryogenesis' mechanism of Cohen, Kaplan, and Nelson. We have assumed that inside 

^Incidentally, note that this is not a general property due to the insertion of QCD sphalerons into the set 
of processes in equilibrium, but is due to the fact that only top Yukawa interactions are fast. If, for instance, 
bottom quark Yukawa interactions were also fast, then we would find that {B + LJ^q is not simply proportional 
to {B — L)' , but depends also on the values of the other charges in equilibrium. 
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the walls of the bubbles nucleated during the electroweak phase transition a space time depen- 
dent 'charge potential' for (partial) fermionic hypercharge is generated. We stress again that 
no discussion about the origin of the charge potential has been given; in particular, since in 
the limit in which all the Yukawa couplings go to zero there is no communication between the 
Higgs and the fermion sectors, we expect that in this limit also the charge potential should go 
to zero. In the traditional approach of CKN there is no trace of this behaviour, and moreover 
we have shown that the results change in a sensible way according to the precise form of the 
charge potential which is considered. We reserve a discussion on this subject for a forthcoming 
publication. Our main interest here was to set a scheme for calculations in the the adiabatic 
scenario in the case in which such a charge potential is present, in order to determine the 
variations in the thermal averages induced by transport effects and the production of baryon 
number by sphalerons inside the bubble walls. 

The main physical point of the paper is that the system should be regarded as a collection 
of subsystems in local equilibrium, in which the thermal averages for the conserved charges are 
not zero but are driven to non vanishing values by transport phenomena. In particular, the 
local equilibrium configurations will correspond to non zero values for {B + L)'. We have deter- 
mined the local equilibrium configuration by means of a chemical potential analysis, calculating 
the values of the thermal averages for the conserved charges by using linear response theory. 
We have considered only the dominant contribution to these averages, in particular, we have 
neglected any fermion mass effect and also the coupling of the Higgs field to the Chern- Simons 
number. 



We find that, in contrast with previous claims fT^, the presence of transport phenomena does 
not prevent baryon number generation inside the bubble walls. The main source of disagreement 
with JPT is the following. In their paper, JPT consider the rate equation in the form 

B = {3fitL + + /^Ti + fJ^i^r), (55) 

where the term on the right hand side has been obtained by considering the variation of the free 
energy of the system due to a 'sphaleron-like' transition involving only the third generation, 
i.e. due to the processes titihiTL ^ and t^bLbii^r ^ 0. Then these authors impose that 
the chemical potential of any particle is proportional to the value of its hypercharge, and so 
they find that the right hand side vanishes as a consequence of the conservation of hypercharge 
(and of fermion hypercharge) in any sphaleronic transition. The point is that, assuming local 
equilibrium of the fast interactions, the chemical potentials of the single particle species are not 
proportional to their hypercharge. In fact, since the single particle numbers are not conserved 
quantities of the system, their abundances are reprocessed by fast interactions as to obtain their 
local equilibrium values. On the other hand, imposing that the particle chemical potentials are 
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proportional to hypercharge, would be equivalent to freeze out any interaction inside the bubble 
wall, both the fast and the slow ones, leaving transport phenomena as the only relevant process. 

Also, if a charge potential for fermionic or total hypercharge is present, transport phenomena 
allow the generation of the baryon asymmetry even in the limit in which the Higgs VEV's go to 
zero. The reason is again that the thermal averages for Q', {B — L)' and BP' are non vanishing 
and then a {B + L) asymmetry can be generated even if the electroweak symmetry is unbroken. 
This is strictly analogous to the well known result of the survival of a B + L asymmetry when 



a. B — L density, eventually of GUT origin, is present |20|. Of course, in the limit in which 



the VEV's go to zero, also the charge potential should go to zero, since no complex phase can 
emerge from the Higgs sector in this case. Then, also this suppression, like the one due to 
vanishing htop, should be made evident by an accurate discussion on the origin of the charge 
potential. As a consequence, our results for the baryon asymmetry, eq. ( ^31) should be probably 
multiplied by a further suppression factor roughly of order {mtop(T)/Ty. 
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